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Symmetric Orbifolds and Entanglement Entropy for Primary Excitations 

in Two Dimensional CFT 
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We use the techniques in symmetric orbifolding to calculate the Entanglement Entropy of a single 
interval in a two dimensional conformal field theory on a circle which is excited to a pure highest 
weight state. This is achieved by calculating the Reney Entropy which is found in terms of a 
2n-point function of primary operators, n being the replica number. 



Entanglement Entropy (EE) has been the subject of 
extensive research in the last few years. Early interests 
in the subject, 000, came through the similarities of 
this quantity with the entropy of black holes through the 
area law [4] [5]. However, EE proved to be a powerful 
tool on its own for understanding the quantum nature 
of physical systems (see for example [6.4S| for reviews and 
references). Apart from being an important quantity in 
quantum information, EE has been introduced as a suit- 
able order parameter for quantum phase transitions. 

When there is entanglement between the degrees of 
freedom of two disjoint regions of space, even if they are 
far apart, measurements in one area affect those in the 
other instantaneously, that is, the effect is not propagated 
by any messenger, rather, it is the result of the quantum 
structure of the system. If an observer confines himself to 
one of the two regions, although he has no access to the 
other, his measurements are affected by that. He makes 
observations in his measurements which are induced by 
the degrees of freedom which are not accessible to him. 
EE is a quantification of his lack of information about 
the subsystem that is accessible to him. 

Lack of information can also be caused by statistical 
distribution of states in a system such as thermal en- 
sembles. In such situations EE will no longer be a use- 
ful measure of quantum entanglement and thus one usu- 
ally studies this quantity when the system is in a pure 
state. Most of the research on this subject has focused 
on the case where this pure state is the ground state of 
the theory. In this article we are interested in excited 
pure states. 

After the proposal of (see also [lol- 12\ for reviews 
and references), which gives a holographic interpretation 
of EE, there has been an even greater interest in the 
subject. The present work grew out of an attempt to 
find the holographic description of the EE of a single 
interval in a two dimensional CFT which has been 
excited by primary operators. The field theory side 
of this problem has already been addressed by two 
different methods in [l3[ and Here we present a 



third method, symmetric orbifolding 15), to address this 
problem which proves useful in finding its holographic 
description [16j . But first some preliminaries. 

Entanglement Entropy in QFT 

Suppose a physical system consists of two subsystems 
A and B and that the whole system is in a pure quantum 
state with the density operator p. Once we take the trace 
of the density operator over the B degrees of freedom, the 
resultant operator is called the reduced density operator 
and is denoted by pa(= trBp). Generically pa will no 
longer be pure and one can associate entropy to it. The 
EE for the subsystem A is denoted by Sa and is de- 
fined as the Von Neuman entropy of the reduced density 
operator, Sa = — tTAPA In Pa- 

A useful mathematical quantity, called the 
Reney Entropy (RE), is defined by the replica trick 



as S A n ^ = j^— \ntrAPA n such that Sa = lina n ^.i S A . 
Here n is a positive integer and is called the replica 
number. In Quantum Field Theory (QFT) this quantity 
can be represented in terms of a path integral. This is 
achieved by considering n copies of the world volume 
of the original theory, A4, and glueing them along the 
entangling subspaces in a cyclic order. This results 
in a space, which we denote by lZ n , and which has 
singularities on the boundaries of entangling subspaces. 
The path integral on 1Z n is denoted by Zn n and is 
defined as 
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This expression is proportional to RE. Except for some 
rare examples it is extremely difficult, if not impossi- 
ble, to calculate Zn n directly. One can go around this 
by transferring the geometric complexities of the world 
volume into the geometry of target space. That is, 
one considers the original nonsingular world volume but 
instead introduces n copies of the target space fields, 
(pi (i = 1,2, ...,n), on that. Instead of gluing the world 
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volumes one now restricts the fields to satisfy certain con- 
ditions along the entangling subspaces 
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where the subscript res stands for restrictions on fields. 
Note that these restrictions replace the nontrivial geome- 
try of lZ n . One way to impose the restrictions is to insert 
the so called twist operators at the boundaries of entan- 
gling subspaces and calculate an unrestricted integral 
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where crj. are the twist operators that enforce restrictions 
through their Operator Product Expansion (OPE) with 
fields. Another way of imposing the restriction is to move 
over to the covering space of the fields, denoted by A4c, 
with a suitable coordinate transformation and perform 
the calculations on this smooth manifold 
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On the covering space the restrictions on fields in the 
integration are taken care of by the geometry of Aic- In 
the following we use this last way of imposing restrictions 
by orbifolding techniques. These restrictions amount 
to identifications of target space fields under subgroups 
of the symmetric group, S n , and hence the name 
symmetric orbifolding. This process has been worked 
out in full detail in [151 ] for two dimensional theories in 
their ground state. The new ingredients in our case are 
the primary operators which excite the theory out of its 
ground state. In what follows, we focus on a generic two 
dimensional conformal field theory and consider a single 
entangling subspace. First a short outline of the method. 

Symmetric Orbifolding 

The outline is as follows [Tij: suppose we start with a 
theory on sphere, parametrized by {z,z), with a flat met- 
ric and with two branch points of order n at u and v. Here 
n is the replica number and u and v are the endpoints of 
the entangling interval. By a coordinate transformation 
to (w(z),w(z)), which behaves as w ~ z x l n at branch 
points, one moves over to the covering sphere with the 
same line element (but a different metric). By a Weyl 
transformation with a conformal factor |^-| 2 , one ends 
up with a third sphere with a fiducial metric ds 2 which 
we have chosen to be flat. 

As the first two spheres are related by diffeomorphism, 
partition functions on the two are equal. This in turn is 
related to the partition function on the third sphere by 
the exponent of the Liouville action imposed by Weyl 
anomaly. A careful calculation of this action results in 
the known expression for RE. 



The case for excitations 

We now wish to calculate the RE for a CFT on a sphere 
with a single branch cut. We further assume that the 
state we start with is a highest weight state with weights 
(h, h). To proceed we use the following parameterisation 
for the sphere 
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We call this the z-sphere for which we have chosen a flat 
metric with a regularisation parameter S. Without loss 
of generality we choose the branch points as following 
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In order for the theory to be in a highest weight exci- 
tation, we create the asymptotic in and out states by 
putting the corresponding operator 0{z, z) at z — z = 
and 0{z,z) at z = z — where 



d(Sj) = 0( Zl z)z 2h z 2jl 5 2 <- h+ V , 



(3) 



and O is a primary operator with weights (h,h). Note 
that the operator O is defined on the north pole cap 
which is parameterised by z and z. One can equivalently 
create the out state at the north pole by introducing the 
adjoint operator as 1 



lim OUz,z) 

z,z— >0 
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lim 0(z,z)z 2h z 2jl S 2 ^ . (4) 



We should now make an n-sheeted Riemann sphere by 
appropriately gluing n copies of the above spheres along 
the branch cuts and calculate the path integral on the re- 
sulting manifold. As stated before, one can equivalently 
calculate a restricted path integral of n copies of the the- 
ory on a single sphere. We denote such a quantity by 
trp^iff) which is defined as 



trpl(6) 
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1 In the following, as a shorthand, we will be rather sloppy with 
this notation and denote the out state by O in the z and z 
argument. So, for example, an out state in the north pole will 
be denoted by O(co). 
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where 9 determines the entangling interval and the ex- 
pression in the denominator is a normalization factor 
which ensures trp® = 1. 

As explained in [HI, the above path integral can be 
properly defined by cutting out circular holes of radius, 
say, e around the branch points as well as the infinity. 
One should then specify the proper boundary conditions 
for fields along the edges of the holes. The precise way 
of performing this procedure has been carried out in [15| 
with full details which we mostly skip. There are two new 
ingredients in our case which should be addressed. The 
first one is that unlike the case of 15[ where the unity 
operator is inserted in z = oo, we insert the primary 
operator O. Secondly, we are also inserting O at z = 0. 
We address these issues in the following. 

Let us define the covering space, the w-sphere, through 
the map 



(G) 



1 _ ( w — 1 \n 
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Near the branch points, the map behaves as z « t n . On 
this sphere the metric is induced through the map as 



ds 2 — — — —dwdw 
dw dw 



(7) 



and the regularisation parameter is also found as 8' — 
a<5sin (0/2) /n. There are several holes on this sphere. 
Two of these are the images of the branch points. There 
are also n holes coming from the images of z = oo and n 
holes from z — 0. 

The prescription of filling these holes should be such 
that upon path integration inside the holes we should 
end up with the desired states at the edges. As for the 
branch points, one such prescription is given in (lBj . This 
is roughly filling the holes with disks of a fiat metric 
which is continuously matched with the metric outside 
the hole over the edge. There is however a curvature 
concentration along the edge which is expected because 
around each branch point there is a cone with an excess 
angle. 

The holes coming from z = oo are naturally filled in 
by the images of \z\ < 1/8 and those coming from z = 
by the images of \z\ < e. The vertex operators O and O 
guarantee the desired wave functional at the edges. 

We now have a closed surface, the w-sphere, on which 
we want to calculate a certain path integral in presence 
of operator insertions. As usual, we choose a fiducial 
metric, which for simplicity, we take it to be a flat metric 
with an arbitrary regularisation parameter 2 , S. This is 
achieved by performing a Weyl transformation on the w- 
sphere, with a factor \dz/dw\ 2 , to make the metric flat. 



We then introduce the coordinate t = w but choose 8 as 
the regularisation parameter which defines t — l/(8 2 t). 

The t-sphcrc thus defined is where we perform our cal- 
culations. Note that the path integrals on the different 
spheres are schematically related as 



[d n tp{z,z)]e- s ^<"< v ^... = J [d V {w,w)}e' s[v] ... 

= e Si f[d^(t,t)}e- s ^... ,(8) 



where dots stand for possible insertions, in case of which, 
the appropriate transformation factors should be in- 
cluded. Sl in the above relation is the Liouville action 
coming from Weyl anomaly and is defined as 



St 



c 
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where c is the central charge of the theory and 



, dz , 



dt 



(9) 



(10) 



The transformation^ turns the restricted path integral 
into an unrestricted one. This is so because the w-sphere, 
being an n-fold cover of the z-sphere, automatically en- 
forces the desired identifications of fields along the branch 
cut. 

As for the insertions, we note that the sequence of the 
Diff x Weyl transformations, keeping the metric invari- 
ant, is a conformal transformation on the whole. There- 
fore for the primary fields O, we have 



and a similar one for O. 

Putting everything together we find that 
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(12) 



where T is the transformation factor for operators, a 
product of those appearing in (JTTJ), and and t' k are 
the images of z = and z = respectively. Z stands for 
partition function and the subscripts z and t denote on 
which sphere the corresponding quantities are calculated. 

We now have all the ingredients to perform our calcu- 
lations. First start with the Liouville part. It turns out 
that there are three different contributions to Sl (see [HI 
for details). One is coming from the kinetic term in the 
1/(5' < \t\ < 1/8 region and the other comes from the 
curvature ring at \t\ = 1/8 3 . These two sum up to an 



We take 8 < S' as a convenient choice for reasons to follow. 



3 In [TJ| these are denoted by S^ 2 ' and S L 3 ^ respectively. 
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expression which only depends on regulators S and d and 
which cancels out in the end. 

The main contribution to Sx, which we call S^, comes 
from the region which is bounded between the edges of 
the holes, on the one hand, and \t\ — 1/5' on the other. 
Given that the metric is flat in this region, the only con- 
tribution to © comes from the Kinetic term which we 
can turn into a boundary term 



bL ~96^ 



dt 



' + c.c. 



(13) 



where the integration is along the boundaries of the re- 
gion. The calculation of this integral for the branch 
points as well as the images of z = oo is identical to 
those in [l5| and nothing changes. As for the images of 
z = 0, it is straightforward to show that near such points 

dz 1 
<t> = log | -77 1 2 « log (c + a(t- t )) , d t <p~ . rr , 
at c + a[t — to) 

where a and c are constants and to is any of the images 
of z — 0. The integral (fl"3]l for these values in the limit 
t — > to will obviously be zero and there is no contribution 
to Sl from the images of z — 0. The upshot is that as 
far as Sl is concerned, operator insertions have no effect. 
Recalling that trp n = e SL Z t /Z™, and consulting (1T2"1) . 
this statement leads to 



trp n (9) 



T 



(UVoO(t k )Q(t' k ))t 
(O(O)O(oo)}" 
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We now find the effect of vertex operators, i.e., cal- 
culate the factor T appearing in (JT2J) . On i-sphere the 
branch points z = (u 7 v) are mapped to t — (— 1, 1). The 
point z = 0, on the other hand, is mapped to n points 
which we denote by tk 



tk 



-i cot 



6 + 2irk 
2n 



fc = 0,l,...,n-l. (15) 



The images of z = oo are denoted by t' k where 



t' k = -icot (— ) , k = 0,l,...,n- 1. (16) 



Note that t' = oo. 

To find T, we need to calculated dt/dz at t = (k — 
0, 1, ...,n— 1), dt/dz &tt = t' k (k — 1, n— 1) and finally 
dt/dz at t = t' . One can then write T as a product of 



n-l 



n-1 



*o fe=l *fe k=0 



(17) 



times a similar expression but with (z, t, h) replaced with 



(z, t, h). One finds that 



dz 
~dl 

dz 
~dl 

dz 
di 



an 



sin(0/2) 



6> + 27rfc 



5 2 asm(9/2) 

k 

,8s2 n 



5' asin(0/2) 



, fc = 0, 1, n — 1 , 
k = 1,2, ...,n- 1 , 

(18) 



which gives the final result as 
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2^) sin(- 1 )(0/2) \ 

(n^oOfa)Q(^)) t 
(O(O)d(oo))" 



2(h+h) 



(19) 



This is our main result. Let us find an approximation to 
this formula in the limit 6 <C 2tt. Recall that 



0(t,*)d(o,o) 



i 



ArA 



t 2h t 2h 



[1 + Qa,a* A < 
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where dots in the second line stand for higher powers of t 
and t and Q A ^ is the operator with smallest dimensions, 
(A, A), in the OPE. This gives 



[[0(t k )0(t k ) = ^^ry sin n {d/2n) ) 
x (1 + O(0( A+A ))) . 



2(h+h) 



We also recall that 

(O(0)6(oo)) a = <5 2( ' i+;i) . 
Putting everything together 

4")(0) = 1 + *±*(I - n)(^) 2 + 0(9^) 



(21) 
(22) 

(23) 



which is in complete agreement with the result in [13j in 
the same limit. 



To further compare our result to those in [13j and [14| , 
we define yet another coordinate s by 



t + 1 
t-1 



(24) 



which maps the t-sphere into a cylinder with unit radius. 
Upon this transformation, the expression in (|19p will find 
the following form 



■4" } (0) 



(0{ e -)6{Q)Y cy 
n- 2 ^ h+h HU n k Zo 0( g± ^)0(^))c 



(25) 
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where in the second line we have scaled the arguments in 
the denominator by a factor n. This is the result found 



[13| and (14 1 . 



EE is now immediately obtained by calculating S = 
d/d n trpb(0) at n = 1. Noting that tr Po (6) = F^\ff) = 
1, we find 



So(0) 



d_ 



s GS {e)~i r ^\e)\ n=l . 



(26) 



13 [3 [II, h is a fair 



Given the Gauge/Gravity duality, 
question to ask for a holographic analogue of the above 
calculations and results. This will amount to identifying 
the gravitational counterparts of the procedure used here. 

In the original proposal of Q, the distinct role of the 
branch points on the boundary result in a desirable 
geometric realisation of EE in terms of minimal lengths. 
This is intuitively understood by extending the curvature 
concentration at the branch points into the bulk. Using 
the covering space, as is the case in this letter, this 
geometric intuition is lost because the singularities at 
branch points are smoothed out. One should then find 
the bulk geometry which corresponds to the covering 
space in field theory. Exciting the theory into primary 
states will then be a matter of turning on appropriate 
bulk fields. This is the subject of flrj ]. 
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